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Abstract 

Background: Social insects, such as honey bees, use molecular, physiological and behavioral responses to combat 
pathogens and parasites. The honey bee genome contains all of the canonical insect immune response pathways, 
and several studies have demonstrated that pathogens can activate expression of immune effectors. Honey bees 
also use behavioral responses, termed social immunity, to collectively defend their hives from pathogens and 
parasites. These responses include hygienic behavior (where workers remove diseased brood) and allo-grooming 
(where workers remove ectoparasites from nestmates). We have previously demonstrated that immunostimulation 
causes changes in the cuticular hydrocarbon profiles of workers, which results in altered worker-worker social 
interactions. Thus, cuticular hydrocarbons may enable workers to identify sick nestmates, and adjust their behavior 
in response. Here, we test the specificity of behavioral, chemical and genomic responses to immunostimulation by 
challenging workers with a panel of different immune stimulants (saline, Sephadex beads and Gram-negative 
bacteria E. coli). 

Results: While only bacteria-injected bees elicited altered behavioral responses from healthy nestmates compared 
to controls, all treatments resulted in significant changes in cuticular hydrocarbon profiles. Immunostimulation 
caused significant changes in expression of hundreds of genes, the majority of which have not been identified as 
members of the canonical immune response pathways. Furthermore, several new candidate genes that may play a 
role in cuticular hydrocarbon biosynthesis were identified. Effects of immune challenge expression of several genes 
involved in immune response, cuticular hydrocarbon biosynthesis, and the Notch signaling pathway were 
confirmed using quantitative real-time PCR. Finally, we identified common genes regulated by pathogen challenge 
in honey bees and other insects. 

Conclusions: These results demonstrate that honey bee genomic responses to immunostimulation are substantially 
broader than the previously identified canonical immune response pathways, and may mediate the behavioral 
changes associated with social immunity by orchestrating changes in chemical signaling. These studies lay the 
groundwork for future research into the genomic responses of honey bees to native honey bee parasites and 
pathogens. 

Keywords: Honey bees, Immunity, Genomics, Social behavior, Social immunity, Chemical ecology, Cuticular 
hydrocarbons, Nestmate recognition 
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Background 

Honey bees are an outstanding model system for study- 
ing the molecular, physiological and social basis of dis- 
ease transmission and resistance. Honey bees are 
plagued by a number of parasites and pathogens 
(reviewed in [1]), and the social colony environment 
(with up to 50,000 densely packed worker bees [2]) pro- 
vides excellent conditions for disease transmission. 
While the innate immune systems of invertebrates and 
vertebrates are surprisingly conserved [3], there can be 
large differences in the numbers of genes involved in the 
different molecular arms of the immune response system 
across species. As is the case with many insects with re- 
cently sequenced genomes [4-8], bees have a much smal- 
ler number of known canonical immune response genes 
relative to Drosophila [9], which is one of the best char- 
acterized models of insect immunity [10]. Thus, bees 
(and other insects) may utilize alternative genetic and 
physiological mechanisms to respond to infections. Fur- 
thermore, honey bees can resist pathogens and parasites 
by employing sophisticated behavioral defense mechan- 
isms, termed "social immunity" (reviewed in [11]). Here, 
we examine the effects of a panel of general immune eli- 
citors (injection with saline, Sephadex beads or bacteria) 
on worker-worker social interactions and chemical com- 
munication, and use whole-genome microarrays to 
characterize global gene expression responses to these 
immune elicitors. 

Honey bee populations have been in decline world- 
wide, with beekeepers recently reporting massive annual 
losses (>30%) (reviewed in [1]). This decline is undoubt- 
edly due in part to the multitude of parasites and patho- 
gens that target honey bees, several of which have only 
recently been identified. Honey bees are host to over 20 
viruses [12-14], as well as a number of bacterial and fun- 
gal pathogens [15-18] including the gut microspordian 
parasites Nosema apis and Nosema ceranae. Honey bees 
are also severely impacted by Varroa mites (Varroa de- 
structor), tracheal mites {Acarapis woodi) and other 
ectoparasites [19]). 

Sequencing of the honey bee genome identified 177 
genes associated with the canonical immune response 
pathways in insects [9]. Innate immune pathways in 
insects, mostly obtained from studies in Drosophila, 
consists of both cellular and humoral responses, which 
can be systemic or local (reviewed in [10,20-24]). Cellu- 
lar immune responses involve a number of differentiated 
hemocytes. Pathogens activate phenoloxidase and asso- 
ciated immune cascades, resulting in phagocytosis, en- 
capsulation, and/or melanization of invading organisms 
or wounds. Humoral responses include cytotoxic mole- 
cules (such as reactive oxygen and nitrogen species), 
lysozymes, cytokines, and antimicrobial peptides 
(AMPs). AMPs are primarily produced by the fat bodies 



in insects. A number of signal transduction pathways are 
involved in moderating immune responses, including the 
JAK/STAT (which seems to respond primarily to tissue 
damage), JNK (which mediates wound repair), Imd 
(which primarily regulates responses to Gram-negative 
bacteria), and Spaetzle/Toll (which generally regulates 
responses to Gram-positive bacteria and fungi) path- 
ways. Activation of one of these signal transduction 
pathways may lead to up- or down- regulation of the 
others depending on the host-parasite system. For ex- 
ample, studies in Drosophila have shown that the rela- 
tionships between immune pathways and effectors can 
be highly pathogen specific and much still remains to be 
characterized [25]. Furthermore, canonical pathways 
derived from Drosophila experiments, may not be fully 
generalizable to other systems. 

Honey bees also employ a number of behavioral 
mechanisms - termed "social immunity" - to reduce the 
impacts of parasites and diseases (reviewed in [11,26]). 
Grooming can remove ectoparasites, and hygienic be- 
havior can reduce levels of brood diseases and depress 
Varroa mite populations (mites feed on developing 
pupae). Workers also coat the interior of their colonies 
with propolis (derived from plant sap), which has anti- 
microbial properties and results in lowered levels of 
immune response genes in individual bees [27]. Further- 
more, diseased or otherwise stressed bees undergo accel- 
erated behavioral maturation, moving from in-hive tasks 
to foraging, which removes them from the brood nest 
and colony, thereby potentially reducing transmission 
[28,29]. Finally, studies in termites suggest that healthy 
nestmates can become "socially vaccinated" by exposure 
to infected nestmates, thereby improving their resistance 
to a pathogen [30]. Future investigations may uncover 
similar defense mechanisms in honey bees. 

Honey bees can also use cuticular hydrocarbons as 
chemical cues to distinguish between healthy and 
immuno-stimulated adult nestmates, and respond differ- 
ently to them [31]. Cuticular hydrocarbons are synthesized 
in oenocytes, which are embedded in the fat body tissue 
under the epithelia, and deposited on the cuticular surface 
[32]. Cuticular hydrocarbon patterns can be modified by 
genotype, physiological state, and environmental context 
(reviewed in [33,34]), including social status in honey bees 
[35]. Immunostimulation of honey bee workers with lipo- 
polysaccharides derived from bacterial cell walls caused 
significant changes in cuticular hydrocarbon profiles in 
worker bees after four hours, and resulted in altered social 
interactions [31]. Healthy bees were more aggressive to- 
wards their nestmates that were coated with extracts of 
immunostimulated nestmates than to nestmates that were 
coated with extracts of healthy nestmates, indicating that 
these changes in chemical profile could alter worker- 
worker interactions. Other studies have demonstrated that 
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Varroa parasitized pupae and adults have modified cuticu- 
lar hydrocarbon profiles [36], while infection with a virus 
can elicit aggressive reactions from healthy nestmates in a 
colony [37]. 

Here, we examined the behavioral, chemical and mo- 
lecular responses in honey bee workers to injection with 
saline, Sephadex beads, and Gram-negative bacteria 
(freeze-killed E. coli cells), six hours after injection. We 
determined if these different immune elicitors stimulate 
unique responses at these three levels; for example, the- 
oretically only E. coli injection should activate the Imd 
pathway, which could lead to distinct chemical profiles 
and behavioral responses. We also determined if immu- 
nostimulation resulted in significantly altered expression 
of previously annotated honey bee immune genes [9], if 
there was overlap with other studies of the effects of 
immunochallenge or parasitization in honey bees and 
Drosophila [38-41], and if there were changes in expres- 
sion of genes associated with cuticular hydrocarbon bio- 
synthesis pathways. 

Results 

Pilot study on the effects of immune elicitors on 
locomotion and mortality 

Lethality of our treatments was assessed 8 hours after 
treatment, in a pilot study using bees from Colony 3. 
After treatment, 1 out of 10 of the control (handled and 
C0 2 anesthetized) and 1 of the 10 bead- or bacteria- 
injected individuals died, respectively. None of the 
saline -injected workers died. 

The behavioral effects of treatments (N = 10 for each 
treatment) were individually assessed on isolated bees 
placed in circular arenas for 5 minutes as in [31]. Statis- 
tical analyses revealed no significant effects of treatments 
on general activity (total time spent self-grooming: con- 
trol: 63 ± 23 sec; saline-injected: 56 ± 15 sec; bead- 
injected: 86 ± 24 sec; bacteria-injected: 83 ± 31 sec; 
Kruskal-Wallis: H(3)=0.58; p=0.8) or locomotion (number 
of lines crossed: control: 50 ± 16 sec; saline-injected: 
29 ± 8.5 sec; bead-injected: 39 ± 12 sec; bacteria-injected: 
43 ± 15.5 sec; Kruskal-Wallis: H(3)=0.507; p = 0.9). 

Effects of immune elicitors on social interactions 

To monitor social interactions, we used a nestmate rec- 
ognition assay in which an individual, 6 hours after 
treatment, was returned to her cage with her original 
nestmates and social interactions were monitored. 975 
bees from Colony 1, separated in 65 petri dishes, were 
used for this experiment. This assay has been used regu- 
larly in bees and ants to examine nestmate recognition 
[31,35]. As in the pilot study with isolated individual 
bees, there was no effect of treatment on locomotion 
(number of lines crossed: control: 6 ± 4.75; saline- 



injected: 4 ± 4; bead-injected: 5.5 ± 3; bacteria-injected: 
2.5 ± 2.5; Kruskal-Wallis: H(2)=-4.6; p>0.05). 

To assess social interactions, agonistic and non- 
agonistic contacts were analyzed by calculating a global 
aggression index (see Table 1). The total aggression 
index was significantly different among the four groups 
(Kruskal-Wallis: H(3)=29.3; p<0.001). Pairwise compari- 
sons revealed a significantly higher aggression index to- 
ward bacteria-injected workers compared to control, 
saline-injected and bead-injected workers (Figure 1A). 

Furthermore, unlike the pilot study with isolated bees, 
the frequency of self-grooming behavior was significantly 
lower in bacteria-injected workers compared to control 
and saline-injected workers, while self-grooming by 
bead-injected workers was intermediate (Kruskal-Wallis: 
H (3, N= 65) = 15,82 p =0.001 followed by post-hoc 
pairwise comparisons p<0.01, Figure IB). However, 
bacteria-injected workers were the targets of significantly 
higher allo-grooming behavior than control workers, 
while saline- and bead-injected workers were intermedi- 
ate (Kruskal-Wallis: H (3, N= 65) =7,97 p =0.046 fol- 
lowed by pairwise comparisons p<0.05, Figure IB). 

Effects of immune elicitors on cuticular hydrocarbon 
profiles 

Non-polar compounds were extracted from honey bee 
workers using pentane, and chemical profiles of cuticular 
extracts were analyzed using gas chromatography (GC). 
In our experimental conditions, all the chemical com- 
pounds we identified by mass spectrometry (MS) ana- 
lysis were hydrocarbons. 

The relative proportions of the cuticular hydrocarbons 
extracted from all four treatment groups from Colony 1 
were significantly different (F (39> 113) =4.9032; p<0.0001, 
Figure 2 A and Table 2). Mahalanobis distances between all 
groups were significantly different (all MD > 8; P < 0.01). 
We obtained similar results for Colony 2 (F (12 , ios) =2.004; 



Table 1 Aggression level index of honey bee behavioral 
interactions 



Abbreviation 


Behaviour Description 


Index 


1 


Immobility 


0 


L 


Locomotion: walking/running 


0 


SG 


Self-grooming 


0 


AG 


Allo-grooming(antennal contact, grooming 
or licking) 


1 


T 


Trophallaxis 


1 


0 


Opening mandibles 


2 


B 


Biting 


3 


F 


Biting with gaster flexion 


4 


S 


Stinging 


5 
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Figure 1 Effects of immunostimulation on behavioral 
interactions. A) Workers are more aggressive towards bacteria- 
injected nestmates. Behavioral responses of worker bees to treated 
nestmates were assayed using a caged nestmate recognition assay 
with bees from Colony 1. Data represents the median ± SIQR 
aggression index. The number of replicates/treatment is 15 control, 
13 saline-, 18 bead-, and 19 bacteria-injected individuals. There were 
significant differences in behavioural responses across the four 
treatment groups (H(3)=29.3; p<0.001). Significant differences in 
post-hoc pairwise comparisons (p<0.05) are denoted by different 
letters. B) Workers engage in more allo-grooming of bacteria- 
injected nestmates. There were significant differences in allo- 
grooming of bees in the four treatment groups (H (3, N= 65) =7.97 
p =0.046). Significant differences in post-hoc pairwise comparisons 
(p<0.05) are denoted by different letters. 

V. / 

p <0.05, Figure 2B and Table 3). Mahalanobis distances be- 
tween all groups were significantly different, except be- 
tween saline- and bead- injected workers (MD<4; p>0.05). 

Discriminant analyses based on the absolute quantities of 
cuticular chemical compounds also revealed significant dif- 
ferences between the four treatment groups (F (51> 102 )=4.32; 
p<0.0001 for Colony 1, Figure 3, and F (45> 89 )=3.03; 
p<0.0001 for Colony 2, data not shown). There did not ap- 
pear to be a significant and consistent effect of treatment 
on the relative proportions or quantities of any of the major 
chemical classes across the two colonies. The total quantity 
of branched alkanes decreased in both colonies in the 
bacteria-injected workers, but this difference was not sig- 
nificant (data not shown). 

Global gene expression responses to immunostimulation 

Microarrays were used to monitor global gene expres- 
sion patterns in the eviscerated abdomens (containing 



epithelial tissue, fat bodies, and oenocytes) of bees from 
the four treatment groups. We examined all transcript ex- 
pression levels across all pairs of treatment groups for sig- 
nificant differences (control x saline, control x bead, 
control x bacteria, saline x bead, saline x bacteria, bead x 
bacteria) and found that 670 unique transcripts were sig- 
nificantly regulated among treatment groups in Colony 1, 
and 1610 unique transcripts were significantly regulated 
among treatment groups in Colony 2 (FDR <0.01, see 
Tables SI and S2 (Additional file 1) for a list of tran- 
scripts). Thus, expression levels for these transcripts were 
significantly different between at least two of the treat- 
ment groups. Of these, 302 transcripts were significantly 
regulated in both colonies (see Additional file 1: Table S3a 
and S3b). This overlap in transcript expression was greater 
than expected by chance (Fishers Exact Test; p<0.001). 
The higher numbers of transcripts in Colony 2 may be 
due to the lower number of biological replicates (4 vs 6 
replicates) or the effect of genetic background, which can 
strongly affect responses to immunostimulation [42,43]. 
Despite some differences in relative expression patterns 
for individual transcripts across colonies, hierarchical clus- 
tering of the 302 common transcripts revealed the same 
overall clustering of treatment groups, (Figure 4). Treated 
groups clustered separately from control groups, and 
saline- and bead-injected groups clustered separately from 
the bacteria- injected group. Relative-fold expression 
values for these transcripts across all of the treatment 
groups can be found in (Additional file 1: Table S3a), and 
the associated p-values for the pairwise comparisons 
between treatments can be found in (Additional file 1: 
Table S3b). 

Effects of individual immune elicitors on gene expression 

Pairwise comparisons identified sets of transcripts differ- 
entially regulated between saline-, bead-, bacteria-injected 
bees and the control bees (FDR < 0.01, Table 4). Overlap 
between colonies for each treatment-control comparison 
was significantly greater than expected by chance (Table 4, 
see Additional file 1: Table S4 for a listing of these tran- 
scripts). Ill, 70, and 117 transcripts, respectively, were 
significantly regulated in these pairwise comparisons in 
both colonies, though these transcripts did not necessarily 
show the same directional patterns of expression between 
colonies. While each treatment resulted in significant ex- 
pression changes in a unique set of transcripts relative to 
controls, there was considerable overlap across treatment 
groups (Figure 5). Indeed, 22 transcripts were significantly 
regulated by all three treatments (Figure 5; see Additional 
file 1: Table S5 for a listing of these transcripts). 

Functional analysis of regulated genes 

Gene ontology analysis of the significantly regulated 302 
transcripts found in both colonies (207 of which had 
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Figure 2 Effects of immunostimulation on the relative proportions of cuticular hydrocarbons. Discriminant analysis was performed on the 
cuticular chemical profiles of control, saline-, bead- and bacteria-injected individuals, based on the relative proportion of each compound. There 
were significant differences across treatment groups in (A) Colony 1, F( 39 _ 113 )=4.9032; p<0.0001), all Mahalanobis distances, p<0.01, and (B) 
Colony 2, F (12 , 10 s)=2.004; p<0.05, all Mahalanobis distances, p<0.05, expect for differences between saline and beads injected individuals, which 
were p=0.4. 



unique Drosophila orthologs with Flybase annotations spermatogenesis, wing disc dorsal/ ventral pattern forma- 
and thus were used in the analysis) revealed an overre- tion, tissue development, post-transcriptional regulation 
presentation of genes involved in immune response, of gene expression and protein polymerization (p<0.05, 
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Table 2 Effects of immunostimulation on cuticular 
hydrocarbon profiles of bees in Colony 1 



JUUjIuI ICC 


Control 


Saline 


Beads 


DCU.ICI la 




Median ± 
SIQR 


Median ± 
SIQR 


Median ± 
SIQR 


Median ± 
SIQR 




Alkanes 






Heneicosane 


U.J I zu.u^ 


n 4Q+n or 


0.55+0.10 


0.43+0.08 


Tricosane 


Z.O 


J.^UIU./ U 


3.70+0.38 


3.11+0.60 


Tetracosane 


U. 1 JXU.UZ 


n 1 7_i_n 1 4 


0.22+0.06 


0.16+0.08 


Pentacosane 


499+0.37 




8.89+4.07 


5.99+2.91 


hexacosane 


u.uuxu.uo 


n R^+n ?r 


0.98+0.19 


0.81+0.26 


Heptac osane 


20.44+2.23 


24.00+4.88 


29.63+4.40 


25.77+8.43 


Octacosane 


0.78±0.05 


0.64+0.10 


0.68+0.06 


0.73+0.08 


Nonacosane 


17.65+1.14 


13.02+2.77 


13.88+3.09 


17.14+2.75 


Triacontane 


0.74±0.07 


0.65+0.21 


0.47+0.17 


0.64+0.15 


Hentricontane 


16.33±2.37 


13.10+4.56 


9.92+4.13 


15.90+3.63 


Tritriacontane 


2.16 


2.03+0.73 


1.35+0.63 


2.17+0.73 


Alkenes 










Tricosene 


0.17±0.03 


0.21+0.07 


0.19+0.06 


0.16+0.04 


Pentacosene 


0.31+0.05 


0.42+0.20 


0.37+0.24 


0.28+0.12 


Heptacosene 




n Qi _i_n C7 


0.72+0.96 


0.51+0.35 


Nonacosene 


1 .08+0.28 


1.41+0.22 


1.23+0.20 


1.16+0.13 


Hentriacontene 
Isomerel 


3.59+0.25 


4.13+0.77 


3.80+0.32 


3.32+0.61 


Hentriacontene 
lsomere2 


5.13+0.30 


5.32+1.11 


4.81+0.54 


4.52+0.63 


Tritriacontene 


12.99+1.19 


13.66+2.58 


1 1.53+1.27 


1 1.1 1+2.13 


Alkynes 










Hentriacontyne 


0.12+0.10 


0.37+0.06 


0.00+0.00 


0.18+0.05 


Tritriacontyne 


1.15+0.28 


1.61+0.32 


1.34+0.32 


1.14+0.21 


Methylalkanes 










11.13- 

Methyl pentacosane 


0.20+0.04 


0.20+0.07 


0.27+0.07 


0.17+0.04 


9.11.13- 

Methylheptacosane 


1.15+0.19 


1 .40+0.25 


1.70+0.40 


1 .20+0.24 


11.13.15- 

Methylnonacosane 


1.06+0.18 


1 .34+0.25 


1.67+0.35 


1.11+0.22 


13.15- 

Methylhentriacontane 


0.62+0.1 1 


0.85+0.15 


1.01+0.21 


0.74+0.14 



Data for this figure were obtained from gas chromatography analysis of 
cuticular washes of worker bees from Colony 1. These data represent the 
relative proportions of each compound found within each of the four 
treatment groups. 



see Additional file 1: Table S6 for a listing of the genes 
in these GO categories). However, none of these categor- 
ies survived the Benjamini correction. Immune genes 
corresponding to the major immune pathways were sig- 
nificantly regulated (defensin-1, relish, domeless, cactus, 
melanization protease 1, death related ced-3/DREDD, 
PGRP-SC2, kayak, spirit among others). Pale {pie), a 
tyrosine hydroxylase involved in melanization and 



wound repair, was also significantly up-regulated by bac- 
teria, bead and saline injection in Colony 1 and bacteria 
and bead injection in Colony 2 [44]. Other genes included 
those involved in cell growth and proliferation (insulin- 
like receptor), cytoskeleton structure (basigin, chickadee, 
twinstar, annexin ix, and isoforms of tubulin), extracellu- 
lar matrix components (pericardin and laminin A), Notch 
signaling {apterous, pebbled, groucho), phagocytosis 
(draper), and cabut, which encodes a transcription factor 
that is regulated by the JNK cascade [45] . 

Of the 22 transcripts (corresponding to 17 unique Fly- 
base genes) significantly regulated by injection with sa- 
line, beads, and bacteria relative to controls, two 
categories were significantly overrepresented: organ 
morphogenesis, p<0.005, and developmental process, 
p<0.001 (see Additional file 1: Table S7 for a listing of 
genes in these GO categories). Again, neither of these 
categories survived the Benjamini correction. Several of 
the previously discussed genes (domeless, insulin-like re- 
ceptor, basigin, twinstar, and groucho) are part of this 
group, as well a serine protease immune response inte- 
grator (spirit) which functions in Toll pathway activation 
[46]. Notably, two genes involved in lipid metabolism 
were also found in this group, which is of particular 
interest since fatty acids are the precursors of cuticular 
hydrocarbon synthesis in insects (reviewed in [47]). 
Bubblegum (bgm) encodes a very long chain fatty acid 
CoA ligase, which plays a role in fatty acid metabolism, 
and was down-regulated in the treatment groups [48]. 
Lipid-storage droplet 2 (Lsd-2) is involved in lipid stor- 
age and accumulation, and was up-regulated in the treat- 
ment groups [49]. Lipid storage droplet 1 (Lsd-1), a 
component of lipid droplets in fat bodies [50], was found 
to be significantly down-regulated in both colonies, but 
not by all treatments. 

GO analysis of the 68 transcripts whose expression was 
specifically altered in bacteria-injected workers (corre- 
sponding to 50 unique Flybase genes) yielded only one sig- 
nificantly overrepresented cluster (immune response; 
p=0.001). These genes included cactus, kayak, draper, 
defensin-1, relish, PGRP-SC2. Other regulated genes 
included those that may be involved in metabolism, such as 
Gr28b, a gustatory receptor that seems to mediate diet- 
related changes in immune response [51], sorbital dehydro- 
genase^ (Sodh-2), an alcohol dehydrogenase of sugars [52], 
and trehalose-6-phospate synthetase 1 (Tpsl) which may 
play a role in protection from hypoxia and anoxic injury 
[53]. See Additional file 1: Tables S8 and S9 for a listing of 
the bacteria-regulated genes and GO categories. 

Comparisons of gene expression patterns with previous 
studies 

We examined overlap in gene expression patterns bet- 
ween our study and previous studies in honey bees and 
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Table 3 Effects of immunostimulation on cuticular 



hydrocarbon profiles of bees in Colony 2 


Substance 


Control 


Saline 


Beads 


Bacteria 




Median ± 
SIQR 


Median ± 
SIQR 


Median ± 
SIQR 


Median ± 
SIQR 


Alkanes 










Heneicosane 


0.84±0.10 


0.98±0.10 


0.82+0.06 


0.87+0.13 


Tricosane 


7.55±3.26 


8.25+1.26 


8.65+1.84 


10.28+2.53 


Tetracosane 


0.45±0.13 


0.46±0.09 


0.46+0.06 


0.54+0.06 


Pentacosane 


1 9.70±4.82 


1 7.89±4.25 


17.94+1.90 


19.85+2.57 


hexacosane 


1.12±0.11 


1 .23+0.1 1 


1 .20+0.04 


1.09+0.14 


heptacosane 


27.03±3.75 


29.08+2.38 


28.24+1.58 


27.46+3.37 


Octacosane 


0.52±0.08 


0.57+0.07 


0.63+0.06 


0.54+0.04 


Nonacosane 


7.82±2.12 


8.52+1.58 


9.11+1.18 


7.71 + 1.52 


Hentricontane 


6.09±1.29 


6.81+1.34 


6.25+1.67 


6.64+1.80 


Tritriacontane 


0.89±0.23 


1.25+0.27 


0.93+0.37 


1.16+0.36 


Alkenes 










Tricosene 


0.93±0.43 


0.79+0.20 


0.76+0.27 


1 .07+0.35 


Pentacosene 


2.30±0.81 


2.28+0.46 


2.24+0.40 


2.53+0.47 


Heptacosene 


0.49±0.07 


0.51+0.12 


0.55+0.17 


0.50+0.13 


Nonacosene 


1 .70±0.26 


1 .94+0.20 


1.81+0.33 


1.55+0.19 


Hentriacontene 
Isomere 2 


2.89±0.78 


3.51+0.86 


2.64+0.83 


2.68+0.42 


Tritriacontene 


7.12+1.27 


7.65+1.43 


5.96+1.04 


6.41 + 1.34 


Alkynes 










Heptacosyne 


1 .93±0.23 


2.18+0.41 


2.14+0.26 


2.08+0.27 


Tritriacontyne 


1 .04±0.23 


1.08+0.19 


1 .00+0.22 


0.80+0.08 


Methylalkanes 










9.11.13- 

Methylheptacosane 


1.53+0.25 


1.59+0.18 


1.48+0.18 


1 .30+0.22 


11.13.15- 

Methylnonacosane 


1.61 ±0.30 


1.70+0.17 


1.56+0.12 


1 .42+0.25 


13.15- 

Methylhentriacontane 


0.92±0.19 


0.97+0.1 1 


0.91+0.07 


0.82+0.15 



Data for this figure were obtained from gas chromatography analysis of 
cuticular washes of worker bees from Colony 2. These data represent the 
relative proportions of each compound found within each of the four 
treatment groups. 



Drosophila. It is important to note that there were large 
differences in study design (including pathogen chal- 
lenge, timecourse and tissue) as well as gene expression 
platform and statistical analyses, which makes direct 
comparisons challenging. However we sought to identify 
common regulated genes that could indicate conserved 
immune response mechanisms across different host 
organisms and pathogen challenges. 

Fourteen genes were found in common between the 
significantly regulated 302 transcripts (239 of which 
were annotated honey bee genes with GB identifiers) 



and the canonical immune response genes annotated 
from the honey bee genome by Evans and colleagues [9] 
(see Additional file 1: Table S10 for a listing of these 
genes). This overlap was significantly greater than 
expected by chance (Fishers Exact Test, p<0.0001; 
Table 5). However, clearly the majority of significantly 
regulated genes are not members of these canonical im- 
mune response pathways. 

Navajas et al. [40] performed a microarray analysis of 
control and Varroa mite-parasitized honey bee pupae 
from two strains (resistant and sensitive to Varroa mite 
infestation). Expression levels of 32 genes (with 22 
matching unique GB identifiers) were significantly 
altered by Varroa parasitization. There was no signifi- 
cant overlap with the 302 transcripts (matched to 
239 GB identifiers) in the current study (Fishers Exact 
Test, p>0.09; Table 5). However, two genes, including 
pie, were regulated in both studies, though directionality 
was not necessarily conserved between studies. 

In a separate study, Alaux and colleagues examined 
the effects of Varroa parasitization and nutrition on 
honey bee worker gene expression using an RNA-seq 
approach [38]. We compared our 302 regulated tran- 
scripts (matched to 239 GB identifiers) with the list of 
genes that were significantly, differentially regulated by 
mite parasitization in pollen-fed bees and found 117 
overlapping genes, including genes involved in immune 
and metabolic processes. Though directionality was not 
necessarily conserved between studies, the overlap was 
greater than expected by chance (Fishers Exact Test, 
p<0.0001; Table 5, see Additional file 1: Table Sll for a 
listing of overlapping genes). Bgm, Lsd-2, and domeless 
were among the genes significantly regulated in both 
studies. A GO analysis of the overlapping transcripts 
(corresponding to 102 unique Flybase genes) found one 
significantly overrepresented category which did not sur- 
vive Benjamini correction (metabolic process, p<0.04). 

Finally, we compared the 302 significantly regulated 
transcripts in our study (matched to 207 unique Flybase 
identifiers) with significantly regulated genes in two sep- 
arate studies examining immune response in fruit flies. 
De Gregario and colleagues [39] identified 400 signifi- 
cantly regulated genes after septic injury of male Dros- 
ophila with bacteria contaminated needles or feeding 
with fungal spores (corresponding to 497 unique Flybase 
genes). Roxstrom-Linquist and colleagues [41] identified 
approximately 390 upregulated genes (corresponding to 
464 unique Flybase genes) in Drosophila melanogaster 
orally infected with protozoa, viruses, bacteria or fungi. 
Eight genes were found in all three studies (Additional 
file 1: Table SI 2), including a number of genes with im- 
mune related functions: pie, cactus, defensin-1, relish 
and PGRP-SC2. Serpin 28D (CG7219J was also regulated 
in all three studies, and is involved in melanization [54]. 
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Figure 3 Effects of immunostimulation on the absolute quantities of cuticular hydrocarbons. Discriminant analysis was used to determine 
if there were significant differences in the cuticular hydrocarbon profiles of control, saline-, bead- and bacteria-injected individuals, based on 
absolute quantity of each compound. There were significant differences across all treatments F(51, 102)=4.32; p<0.0001. 



Quantitative real-time PCR validation of expression of 
candidate genes 

We examined expression levels of several candidate 
genes identified in the microarray study in a third bio- 
logical replicate (Figure 6). Three treatment groups were 
used: control, saline-injected, and bacteria-injected 
worker bees. As in the microarray study, defensin 1 
levels increased with treatment, and were significantly 
higher in bacteria-injected bees relative to controls, and 
intermediate in saline- injected bees (F(2,24)=6.49, 
p=0.0056, see Figure 6 for results of Tukey HSD post- 
hoc pairwise comparisons and Additional file 1: Table S3a 
for expression levels obtained from the microarray ana- 
lysis). As in the microarray study, levels of bubblegum 
decreased significantly with treatment (F(2,24) =32.56, 
p<0.0001), while levels of Lsd-2 (F(2,24)=392.1, p<0.0001) 
and pale increased (F(2,24)= 27.98, p<0.0001). Similar 
to the results of the microarray analysis, the effects of 
treatment on expression of three genes involved in 
Notch signaling were less dramatic, but nonetheless 
expression of apterous (F(2,24)=6.76, p=0.0047), grou- 
cho (F(2,24)=8.44, p=0.0017) and pebbled (F(2,24) 
=4.41, p=0.023) were significantly affected by treat- 
ment. Expression levels were higher in saline-injected 
bees than controls for all three genes, and higher in 
bacteria-injected bees relative to control bees for 
apterous. 



Discussion 

We have demonstrated that immunostimulation of adult 
honey bee workers results in altered social interactions 
only in the case of bacteria-injected bees. However, all 
the treatment groups have significantly different cuticu- 
lar hydrocarbon chemical profiles, and all treatments 
caused large scale changes in gene expression patterns 
in abdominal fat body and epithelial tissue, which 
includes oenocytes. It is important to note that these 
dramatic changes occurred in a relatively short time- 
frame, only six hours after immunostimulation. Import- 
antly, this study demonstrates that expression of a large 
number of genes, not just those canonically associated 
with immune response pathways, are modulated by 
immunostimulation, as has been demonstrated by sev- 
eral other studies of genome-wide responses to immune 
challenges (for example, [38-41]). However, comparisons 
with studies in Drosophila found limited overlap of gene 
expression patterns, and the handful of commonly regu- 
lated genes were key members of canonical immune re- 
sponse pathways. 

It is somewhat perplexing that only bacteria-injected 
bees were subjected to significantly increased allo- 
grooming and aggression. All treatments caused changes 
in chemical profiles, and the changes did not appear to 
be substantially larger in bacteria-injected bees, though 
there was a trend for lower total quantity of branched 
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Figure 4 Hierarchical clustering of significantly regulated 
genes. We performed hierarchical clustering analysis on the 302 
significantly, differentially expressed transcripts (FDR<0.01) similarly 
regulated in Colony 1 (A) and Colony 2 (B). Both colonies 
demonstrated the same overall grouping for experimental 
treatments: saline (S) and bead (B) injected bees formed a sister 
group and were closer in transcript expression to bacteria (Bac) 
injected bees than controls (C). Control and bacteria-injected bees 
had the most disparate transcript expression. Colours denote 
differences in log2 expression relative to the mean expression across 
the four treatment groups, according to the scale shown. 



alkanes in this group. However, previous studies have 
suggested that alkenes, rather than alkanes, play a sig- 
nificant role in nestmate recognition [55]. Based on the 
hierarchical clustering analysis, the gene expression pat- 
terns of bacteria-injected bees were clearly distinct from 
saline- and bead-injected bees. However, overall it did 
not appear that bacterial injection altered expression of 
an entirely unique set of genes (out of 302 transcripts, 
expression of 68 transcripts were significantly regulated 
only by bacterial injection, relative to controls, and ex- 
pression of many of these genes changed non- 
significantly in the other treatments as well). Thus, per- 
haps bacteria injection did not cause change in expres- 
sion of a large, unique subset genes, but rather caused 
more substantial changes in the magnitude of expression 
levels of some genes, as indicated by hierarchical cluster- 
ing. Alternatively, it may be that expression of a few key 
genes and/or chemical profile components in bacteria- 
treated bees provoked behavioral changes in nestmates. 

It remains to be determined if infection with honey 
bee specific parasites and pathogens causes similar 
changes in gene expression and behavioral responses, 
though there is some indication that infection with 
Nosema microsporidia does stimulate expression chan- 
ges in a significantly overlapping set of genes (Holt, 
Aronstein and Grozinger, unpublished data). It also re- 
mains to be determined if these behavioral changes are 
adaptive: it is possible that bees have evolved to "match" 
the strength of the signal to the virulence of the patho- 
gen, and thus infection with bacteria would elicit a 
greater change in the behavioral responses of nestmates 
than simple cuticular wounding. However, colony-level 
assays would need to be performed in order to deter- 
mine if these changes in social interactions actually re- 
sult in differences in the spread of pathogens through 
the colony and their impacts on infected individuals. For 
example, the increased grooming behavior we observed 
towards bacteria-injected individuals could facilitate the 
spread of pathogens through the colony (as observed in 
carpenter ants, [56]), slow the spread by causing isola- 
tion or removal of infected individuals [56,57] or reduce 
the impact on the infected individual [58,59]. 



Richard et al. BMC Genomics 2012, 13:558 
http://www.biomedcentral.com/1471 -21 64/1 3/558 



Page 10 of 17 



Table 4 Analysis of overlap among treatment groups 



Treatment # of significantly, differentially 

group regulated transcripts 





Colony 1 


Colony 2 


In both 


Saline x Control 


340 


555 


111 


Bead x Control 


213 


398 


70 


Bacteria x Control 


278 


1198 


117 


Total 


*557 


*1453 


*206 



p- value 
(Fisher's 
Exact Test) 



p < 0.001 
p < 0.001 
p < 0.001 



Fisher exact tests found more overlap than expected (p<0.001) for 
significantly, differentially expressed transcripts (FDR<0.01) in the saline-, bead- 
and bacteria-injected vs control groups between the two colonies, 
^duplicate transcripts between groups removed. 



Our studies also demonstrate that immunostimulation 
elicits complex gene expression changes in the epithelial 
tissues in honey bee workers. Immunostimulation 
resulted in significant expression changes of 302 com- 
mon transcripts in worker bees from the two colonies 
examined. Only 14 of these genes corresponded to previ- 
ously annotated immune genes identified from the 
honey bee genome [9] . Several other biological processes 
were modified, including cell growth and proliferation, 
cytoskelatal structure, metabolism and components of 
the Notch signaling pathway. Cell growth, proliferation, 
and migration, particularly involving actin-mediated 
cytoskeletal changes, are required for repairing epithelial 
wounds [24]. The Notch signaling pathway has not yet 
been linked to immune response or wound repair, but 
wound repair uses many of the same developmental 
pathways that function during dorsal closure in Drosophila 
development, which does involve Notch signaling [60]. 
Alternatively, since the insect fat body is involved in 
regulating many key processes, including metabolism 
(which was commonly regulated in our study and by 
Varroa parasitization [38]), these changes may reflect 
general physiological changes after immunostimulation 
or stress [61]. Expression changes in three genes of the 
Notch signaling pathway {apterous, groucho, and pebbled) 
were confirmed using quantitative real-time PCR. Inter- 
estingly, expression was significantly higher in saline- 
injected bees relative to controls for all three genes, but 
only expression of apterous was affected in bacteria- 
injected bees, suggesting that changes in Notch signa- 
ling may be modulated temporally or by other signaling 
pathways. 

We found significant expression changes of a number 
of key immune response genes (for a review of the func- 
tion of these genes, see [10,23,24]). The JAK/STAT path- 
way is regulated by the Domeless receptor; domeless 
expression was significantly regulated by immune stimu- 
lation in both genotypes of bees in our study. Activation 
of the IMD pathway requires cleavage of Relish by 
the caspase DREDD; both Dredd and Relish were sig- 
nificantly regulated in our study. We also observed 



Bac vs C 



Be vsC 




Sal vs C 

Figure 5 Effects of specific immunostimulants on gene 
expression. Pairwise comparisons identified sets of transcripts 
differentially regulated between saline-, bead-, bacteria-injected bees 
and control bees in both colonies. A Venn diagram demonstrates 
that there is considerable overlap in the effects of each treatment, 
but also treatment-specific effects on gene expression. The numbers 
represent the number of genes in each category. 



significant regulation of PGRP-SC2 which suppresses ac- 
tivation the IMD pathway, with bacteria-injected bees 
showing the highest levels of PGRP-SC2 expression (data 
not shown). In fruit flies, PGRP-SC1 and PGR-SC2 may 
function in preventing over- activation of the IMD path- 
way [62]. The IMD pathway triggers the JNK pathway, 
which activates the transcriptional regulator AP-1 
(which contains Kayak/D-fos), and AP-1 in turn nega- 
tively regulates Relish-dependent transcription. We 
found significant regulation of kayak in our study. Inter- 
estingly, we also found significant changes in expression 
of cabut, which is regulated by the JNK pathway but has 
not yet been linked to immune function [45]. The Toll 
pathway operates through transduction factors including 
spirit, which acts extracellularly and upstream of 
Spaetzle [46] and the NF-kB protein Dorsal (which is 
negatively regulated by cactus). We found significant 
regulation of spirit and cactus in our study. Pale, which 
plays an important role in melanization and wound re- 
pair, and draper, which functions in phagocytosis, were 
also significantly regulated in our study. As a Gram- 
negative bacteria, it would be expected that E. coli would 
primarily stimulate activation of the IMD pathway. How- 
ever, we observed changes in gene expression of members 
of the Toll pathway, including cactus, spirit and defensin- 
1, which exhibits Gram-positive antimicrobial activity and 
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Table 5 Comparisons of significantly regulated genes 



with previous studies 



Gene list 


# of 
significant 
transcripts 


# of transcripts 

in common 
between studies 


p-value 
(Fisher's 
Exact Test) 


Evans et al. (2006) 


166 


14 


p < 0.0001 


Navajas et al. (2008) 


22 


2 


p > 0.09 


Alaux etal. (2011) 


3704 


117 


p < 0.0001 



The 302 genes that were significantly regulated by immunostimulation in both 
colonies in our study were matched to 239 GB numbers and compared with 
the immune genes identified during annotation of the genome (Evans et al., 
2006), and genes significantly regulated by Varroa parasitization of honey bees 
(Navajas et al., 2008 and Alaux et al., 201 1). A Fischer's Exact Test was used to 
determine if this overlap was significant, assuming that all genes present in 
the background lists for the arrays were in common between studies. In order 
to ascertain overlap, only genes in previous studies that were present on our 
array platform were included. 

is not regulated by the IMD pathway [63]. Thus, there is 
likely considerable cross-talk between the pathways. 

Our studies also identified several genes which may 
play a role in altering cuticular hydrocarbon patterns. 
Cuticular hydrocarbons are synthesized primarily in the 
oenocytes (reviewed in [47]), which are embedded in the 
fat body of adult honey bees [64]. Cuticular hydrocarbon 
biosynthesis involves activation of fatty acids by an acyl- 
CoA synthetase, chain elongations of fatty-acyl-CoAs to 
produce very long chain fatty acids, and subsequent con- 
version to a hydrocarbon, likely by a p450 enzyme [47]. 
Fatty acids are stored in lipid droplets in the adipoctye 
cells of the insect fat body [61]. Fatty acids can be 
released from droplets in the adipocytes and accumulate 
in the oenocytes; this occurs under starvation conditions 
in particular [65]. This process is mediated in part by 
lipid storage droplet-2 (Lsd-2): increased Lsd-2 expres- 
sion in the fat bodies decreases lipid movement to the 
oenocytes. We found increased expression of Lsd-2 in 
immunostimulated bees (see Figure 6), suggesting 
reduced movement of lipids to oenocytes, and perhaps 
reduced levels of cuticular hydrocarbons. We did ob- 
serve a decrease in the total relative quantity of all 
branched alkanes in bacteria-injected workers in both 
colonies but this difference was not significant. 

Bubblegum (bgm) activates long chain fatty acids to 
form acyl-CoAs (reviewed in [66]), a key step in cuticu- 
lar hydrocarbon biosynthesis. Bgm was originally 
described as a Drosophila mutant that resulted in ele- 
vated levels of very long chain fatty acids and neurode- 
generation [48]. Bgm homologs have been identified in 
numerous species, including humans and mice, and have 
been demonstrated to activate long chain (CI 6) and very 
long chain (C24) fatty acids [67]. In our study, bgm ex- 
pression was significantly decreased relative to controls 
(see Figure 6). 

Despite large differences in study designs and analysis 
methods, we found some overlap in gene expression 



with previous studies examining the effects of Varroa 
mite parasitization on honey bees [38,40]. Varroa-re- 
sponsive genes were significantly associated with basic 
cellular processes, including cell organization, biogenesis 
and metabolism. We also found significant changes in 
functional categories associated with basic cellular pro- 
cesses, such as cell growth, proliferation and cytoskeletal 
structure. Varroa parasitization also caused changes in 
expression of pale. As discussed above, this may repre- 
sent cellular mechanisms for wound-healing. Expression 
of potential hydrocarbon synthesis genes, namely bgm 
and Lsd-2, were also regulated by Varroa parasitization. 
Indeed, Varroa parasitized pupae and adults have modi- 
fied cuticular hydrocarbon profiles [36]. These differ- 
ences are likely responsible for stimulating hygienic 
behavior, in which diseased larvae are removed by adult 
worker bees, a key component of Varroa resistance [11]. 

Comparison with two previous studies [39,41] examin- 
ing the effects of immunostimulation on Drosophila glo- 
bal gene expression patterns revealed conserved changes 
in expression of key immune genes in Drosophila and 
honey bees (including relish, cactus, defensin-1, spirit, 
PGRP-SC2, and pale), but otherwise limited overlap in 
the significantly regulated genes. The lack of similarity 
could represent species-specific immune responses or 
simply technical differences - for example, Roxstrom- 
Lindquist [41] orally infected young male flies with bac- 
teria, fungi and microsporidia, and measured whole- 
body gene expression changes in only two replicates 
using Affymetrix microarrays. 

Conclusions 

Our results suggest that immunostimulation of honey bee 
workers causes significant changes in gene expression pat- 
terns, cuticular hydrocarbon chemical profiles, and, in the 
case of bacteria-injection, social behavior. While bacteria- 
injection did not cause expression changes in an entirely 
unique set of genes or production of unique cuticular 
hydrocarbons compared to other treatment groups, it 
does appear that the magnitude of the expression and pro- 
duction changes was greater, and this may have resulted 
in the altered behavioral responses. As demonstrated by 
other studies examining genome-wide expression changes 
associated with immune challenges [38-41], we found that 
even a short-term immune challenge can result in dra- 
matic changes in gene expression that encompass far 
more genes than those represented by canonical immune 
response pathways. For example, we found several genes 
associated with the Notch signaling pathway, which sug- 
gests this pathway may also play a role in mediating im- 
mune responses. Furthermore, our study has highlighted 
potentially new candidate genes for regulating cuticular 
hydrocarbon synthesis in insects. These chemicals serve 
many functions in insects, including operating as sex and 
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Figure 6 Quantitative real-time PCR validation of expression patterns of candidate genes. Expression levels of seven candidate genes 
(relative to actin) were analysed in a third biological replicate using quantitative real-time PCR. Mean expression levels for each treatment group 
are normalized to expression in the control treatment group, for graphical representation. Significant differences in expression levels across the 
three treatment groups were determined using an ANOVA with treatment as a variable, followed by post-hoc Tukey HSD pairwise comparisons. 
Different letters denote significant differences in expression (p<0.05). Nine, eight, and ten individual bees were used for the control, saline- 
injected, and bacteria-injected treatment groups, respectively. 



caste-specific pheromones, and most research on their 
biosynthesis has focused on desaturase enzymes and 
p450s [34]. These studies lay the groundwork for future 
work examining the molecular pathways that mediate im- 
mune responses to acute stimulators and chronic infec- 
tions in bees and other insects, and has paved the way for 
research into the genes that regulate social immunity. 

Methods 

Honey bee stocks 

Honey bee colonies were maintained according to stand- 
ard practices at the Lake Wheeler Honey Bee Research Fa- 
cility at North Carolina State University in 2008. Workers 
for these studies were obtained from three source col- 
onies, headed by queens each instrumentally inseminated 
with semen from a single, different, male (Glenn Apiaries, 
Fallbrook, CA). The specific source colonies used for the 
different experiments were Colonies 1, 2 and 3, and are 
listed in the experimental details below. Since honey bees 
are haplodiploid, the coefficient of relatedness among the 
workers in each colony was 0.75, thereby reducing vari- 
ation in gene expression responses, chemical profiles, and 
presumably behavioral responses among the nestmates. 
Honeycomb frames of late-stage pupae were removed 
from the colonies and placed in incubators overnight at 
33°C, 50% relative humidity (RH). 

Honey bee rearing 

Newly emerged bees (<12 hours old) were brushed from 
the frames and placed into modified 10 cm Petri dishes 



in groups of 15. Dishes were maintained under red light 
in a temperature and humidity-controlled environmental 
room (33°C, 50% RH, Phytotron Facility, NCSU). Bees 
were fed 50% sucrose/water solutions and 45% honey/ 
45% pollen/ 10% water paste ad libitum. Food was 
replaced every two days. Bees were also exposed to 0.1 
queen bee equivalents of queen mandibular pheromone 
(QMP) (Pherotech, Vancouver, Canada). Every day, 10 ul 
of QMP (0.01 queen equivalents/ uL in an isopropanol/ 
1% water solution) was placed on a microscope slide and 
allowed to evaporate before being placed in the cage. 
This amount of QMP mimics a live queen in assays of 
worker behavior and physiology [68,69] and thus should 
help simulate normal rearing conditions. 

Experimental treatment 

When the bees were 10 days old, five individual bees 
were removed from each cage (leaving 10 in the cage). 
Three of these were then subjected to one of four treat- 
ments. The first group of bees was handled and anesthe- 
tized with C0 2 for 1 minute (control treatment). A 
second group of bees was anesthetized and injected with 
8 ul of sterile bee saline (130 mM NaCl, 6 mM KC1, 
4 mM MgCl 2 , 5 mM CaCl 2 , 160 mM sucrose, 25 mM 
glucose, 10 mM 4-(2-hydroxyethyl)-l-piperazineethane- 
sulfonic acid in distilled water, pH 6.7, 500 mOsmol, as 
in [31]). A third group of bees was anesthetized and 
injected with 8 ul of a CM-25 Sephadex beads (Sigma- 
Aldrich, Steinheim, Germany) solution (0.01 g of beads 
mixed with 500 ul sterile bee saline and vortexed, 
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resulting in approximately 110 beads/8 ul). The fourth 
group of bees was anesthetized and injected with 8 uL of 
an E. coli (JM101, Sigma, St Louis, MO) solution sus- 
pended in sterile bee saline (3.8* 10 5 cells/bee, following 
a protocol modified from [70]). The bacteria were grown 
in LB media (Sigma, St Louis, MO), collected and resus- 
pended in the sterile bee saline at the appropriate con- 
centration (4.75*10 4 cells/uL), and then stored at -80°C 
until it was thawed for the injections. Injections were 
performed into the abdominal cavity through tergites 
with a nano-injector equipped with a glass needle 
(Schley Compact Model II Instrument; Honey Bee In- 
semination Services, Davis, CA, US), using a binocular 
microscope (Leica MZ6 stereomicroscope, Leica Micro- 
systems, Buffalo Grove, IL). Treated bees were marked 
with a dot of Testor s paint on their thorax and main- 
tained individually in 10x10x7 cm 3 Plexigas cages in an 
incubator under red light at 33°C, 50%RH for 6 hours, 
with 50% sucrose. After 6 hours, two treated bees per 
cage were collected onto dry ice and stored at -80°C for 
microarray or chemical analysis, while the third was 
reintroduced (individually) to her original Petri dish 
cages for behavioral analysis. 

Behavioral analysis was performed on bees from Col- 
ony 1 (15 control, 13 saline-, 18 bead-, and 19 bacteria- 
injected individuals), while chemical analysis was per- 
formed on bees from Colony 1 (13 control, 13 saline-, 
12 bead- and 16 bacteria-injected individuals) and 
Colony 2 (12 control, 11 saline-, 13 bead- and 12 
bacteria-injected individuals). Microarray analysis was 
performed on 6 individuals/treatment in Colony 1, and 4 
individuals/treatment in Colony 2. Colony 3 was used 
for pilot studies of the deleterious effects of the treat- 
ments (see below). 

Behavioral assessment of the deleterious effects of 
treatment 

Assessment of grooming behavior and locomotion 

To determine whether the treatments induced any dele- 
terious effects, we conducted a preliminary study in 10- 
day-old honey bees (N = 10 per treatment, from Colony 
3). Immediately following the injection, individual work- 
ers were placed into the Petri dish cages, where they 
were later observed. Grooming behavior and locomotion 
were assessed after four hours. To assess locomotor ac- 
tivity, a cross was drawn under the circular arena and 
the number of lines crossed during 5 minutes was used 
as a locomotion index. All studies were conducted blind 
to the treatment group. Mortality was assessed 8 hours 
after treatment and only two workers died. 

Assessment of social interactions 

Six hours after being exposed to one of the four treat- 
ments, individual workers were reintroduced to their 



original Petri dish cage. Assays were conducted blind to 
the treatment group. Each cage was used only once. Be- 
ginning 30 seconds after the reintroduction, the number 
and type of social behaviors of the non-treated bees to- 
ward the treated bee were recorded through a scanning 
method, every 10 seconds for 5 minutes. Aggression 
levels were defined using standard practices, as 
described in Table 1 [31,35]. For each treated bee, the 
aggression score was calculated by multiplying the num- 
ber of non-treated bees observed performing a particular 
behavior during the observation period by the aggression 
level index (Table 1) for that behavior, and summing 
across all interactions observed. 

Hydrocarbon extraction and analysis 

Cuticular hydrocarbons were extracted by submerging each 
frozen bee individually into 1 ml n-pentane with 10 ug of 
hexadecane (Sigma, St Louis, MO, U.S.A.) as an internal 
standard, and agitating gently for 10 min. Whole bees were 
used for all samples. The solvent was reduced to dryness 
with a gentle stream of N 2 , and the sample resuspended in 
100 ul. Hydrocarbons were quantified by gas chromato- 
graphy and identified by gas-chromatography-mass spec- 
trometry (GC-MS). We analyzed 2 ul of solution on a 
HP5890II GC (Agilent, Palo Alto, CA, U.S.A.) equipped 
with a flame-ionization detector and interfaced with a 
HP ChemStation (Rev. A.09.03). Splitless injection was 
made into a 30 m x 0.32 mm xl um HP-5 capillary co- 
lumn operated at 150°C for 2 min, increased at 15°C/min 
to 250°C then at 3°C/min to 300°C and held at this 
temperature for 15 min. The injector and detector were 
held at 300°C and 310°C, respectively. Compound identi- 
fication was achieved by splitless capillary GC-MS using 
a HP6890 GC and a model 5973A MSD with an elec- 
tron impact ion source and a HP-5 capillary column 
(30 m x 0.25 mm ID x 0.25 um film thickness). This 
protocol was modified from [31]. 

Statistical analyses of behavioral and chemical studies 

For behavioral assays, owing to sample size, nonpara- 
metric procedures were used to assess significant va- 
riations in our data. Non-parametric Kruskal-Wallis 
ANOVAs on ranks for global comparison were per- 
formed, followed by corrected Mann-Whitney U tests 
for multiple comparisons to compare significant differ- 
ences in aggression indices between treatment groups. 
For statistical analyses of the chemical profiles, all com- 
pounds present above 5% relative levels were used. To 
assess the similarities of cuticular hydrocarbon profiles, 
a stepwise discriminant analysis of the relative propor- 
tions of the compounds was employed using Statistica 
6.0 (StatSoft, Inc., Tulsa, OK, U.S.A.) as in [71]. Prior to 
analysis, each peak area was standardized according to 
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Table 6 Primer sequences used to quantify expression of 
candidate genes from a third biological replicate using 
quantitative real-time PCR 

Accession number/ 
reference 

XM_003249 187.1 



Gene name 


Primer 
direction 


Sequence 


Lipid Storage 
Droplet 2 


Forward 


AATCGTGCTCG 
CAAAATACC 




Reverse 


TOGCAATAGGT 
AATGAI I I I ICA 


Bubblegum 


Forward 


^GCAAATCT 
TGGTGCAA 




Reverse 


TCTOTGGm 
GCOTCGT 


Pale 


Forward 


TAOTCGTCGCG 
GATAGOT 




Reverse 


TACGCGAACGAT 
GTOTCAG 


Groucho 


Forward 


ACTGTGGACAGG 
AGGATOG 




Reverse 


ATGAAGGACCTC 
GACGmG 


Draper 


Forward 


ATCCTGCCACTG 
GTTATTGC 




Reverse 


AGTCAACTCCCT 
CCCAACCT 


Apterous 


Forward 


AGGAAGAGGAAA 
CCCAAGGA 




Reverse 


TCTGCGAGAGTO 
C1TCAGA 


Pebbled 


Forward 


CGCTGATACGTCA 
CCTCAGA 




Reverse 


GCGAAC1TCTCTC 
CACGTO 


Defensin 1 


Forward 


TGCGCTGCTAACT 
GTCTCAG 




Reverse 


CGAAACGmGTC 
CCAGAG 


Actin 


Forward 


CCTAGCACCATCC 
ACCATGAA 




Reverse 


GAAGCAAGAATO 
ACCCACCAA 



XM_624222.3 

NM_001 01 1633.1 

XM_392425.3 

XM_624852.2. 

XM_392622.4 

XM_003249649.1 

(Evans, 2006) 
NM_00101 1616.2 
(Huising & Flik, 2005) 



Primer sequences were designed in Primer3, v. 4.0 (http://frodo.wi.mit.edu/) or 
obtained from previous publications. 



[72]. Multivariate analysis of variance is known to be re- 
sistant to non-normality in the data because of skewness. 
However, before conducting the analyses, we tested for 
normality and equality of variance over all groups, both 
of which were non-significant. 

Microarrays 

Individual bees were thawed and dissected under cold 
RNAlater (Qiagen, Valencia, CA). Abdomens were evis- 
cerated and the cuticles, with the associated fat bodies 
and oenocytes, were collected. RNA was extracted using 
an RNeasy kit (Qiagen). RNA was quantified using a 



Nanodrop 1000 spectrophotometer (Thermo Fischer 
Scientific, Wilmington, DE) and sample integrity and 
quality was monitored using agarose gel electrophoresis 
to confirm the presence of ribosomal RNA bands. 
500 ng of RNA/individual was amplified using the 
Ambion MessageAmp II aRNA Amplification kit 
(AM1751, Life Technologies, Grand Island, NY). 5 ug of 
amplified RNA from each sample were labeled inde- 
pendently with Cy3 and Cy5 dyes using the ULS aRNA 
fluorescent labeling kit (EA-006, Kreatech, Amesterdam, 
Netherlands). Samples were hybridized to microarrays 
(two samples/array) in a loop design with dye swaps 
incorporated, using 24 microarrays for Colony 1 (6 bio- 
logical replicates, 2 technical replicates per sample) and 
16 arrays for Colony 2 (4 biological replicates, 2 tech- 
nical replicates per sample). Whole genome microarrays 
(<http://www.ebi.ac.uk/arrayexpress/arrays/A-MEXP-755>) 
were purchased from the W.M. Keck Center for Func- 
tional Genomics at the University of Illinois, Urbana- 
Champaign. Arrays were scanned using the Axon Gene- 
pix 4000B scanner (Molecular Devices, Sunnyvale, CA) 
using GENEPIX software (Agilent Technologies, Santa 
Clara, CA). 

Analysis of gene expression 

Any spots with an intensity of less than 100 (the back- 
ground level on the arrays) were removed from the ana- 
lysis. Also, spots present on less than 12 out of 24 arrays 
for the first colony and 8 out of the 16 arrays for the sec- 
ond colony were excluded from the data set as well. Ex- 
pression data was log-transformed and normalized using 
a mixed-model ANOVA (proc MIXED, SAS, Cary, NC) 
with the following model: 

Y = [A + dye + block + array* dye + array* block + e 

where Y is expression, dye and block are fixed effects, 
and array, array* dye and array*block are random effects. 
Transcripts with significant expression differences be- 
tween groups were detected by using a mixed-model 
ANOVA with the model: 

Y = fi + treatment + spot + dye + array + e 

where Y represents the residual from the previous 
model. Treatment, spot and dye are fixed effects and 
array is a random effect, p-values were corrected for 
multiple testing using a false discovery rate < 0.01 (proc 
MULTTEST, SAS). Hierarchical clustering, using the 
ward method was performed in JMP 9.0.2 (SAS, Cary, 
NC). Gene ontology analysis was performed using 
DAVID version 6.7 [73,74] with a cutoff of p <0.05. For 
all gene ontology (GO) analyses, array transcripts were 
matched to their Drosophila orthologs in Flybase (http:// 
flybase.org/). All of the array transcripts with Drosophila 
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orthologs were used as a background list. Analysis of 
overlap between the significantly regulated transcripts in 
our study or between significantly regulated genes in our 
study and other gene expression studies was performed 
using Fisher s Exact Tests (Dr. Oyvind Langsrud, Statistics 
Norway, <http://www.langsrud.com/fisher.htm>). Com- 
mon, differentially expressed transcripts or genes between 
studies were identified with Venny [75] . For overlap com- 
parisons between our study and other studies examining 
gene expression in honey bees ([9,38,40]), we converted 
microarray transcript identifiers (AM numbers) to GB 
identifiers from the honey bee genome annotation (see 
BeeBase http://hymenopteragenome.org/beebase/). For 
comparisons between our study and other studies examin- 
ing gene expression in fruit flies ([39,41]), we converted all 
genes to Flybase identifiers to find common overlap. The 
array data has been deposited on the ArrayExpress website 
(E-MEXP-3708). 

Validation of candidate gene expression patterns using 
quantitative real-time PCR 

Worker bees from a fourth source colony were reared as 
before. Control, saline-injected, and bacteria-injected 
treatment groups were generated (with 9, 8 and 10 indi- 
vidual bees per treatment group, respectively) as before, 
but the injected E. coli bacteria were collected from an 
actively growing culture just prior to injection. Bees were 
collected onto dry ice four hours after treatment and 
stored at -80°C. RNA was extracted from eviscerated 
abdomens of individual bees as above. cDNA synthesis 
and quantitative real-time PCR was performed as in 
[76]. Expression levels of candidate genes were normal- 
ized to actin. Significant differences in expression levels 
among treatment groups were determined using an 
ANOVA followed by post-hoc pairwise comparisons 
with Tukey HSD tests (JMP 9.0.2, SAS, Cary NC). Pri- 
mer sequences are given in Table 6. 

Additional file 



transcripts regulated by all treatments in Colonies 1 and 2 (FDR < 0.01). 
Table S8 The 68 significantly, differently expressed transcripts in bacterial 
x control treatment groups in common between Colonies 1 and 2 
(FDR<0.01). Table S9 Overrepresented GO categories (p < 0.05) amongst 
the 68 transcripts regulated by bacterial treatment only in Colonies 1 and 
2 (FDR < 0.01). Table S10 14 genes that overlap between canonical 
immune genes (Evans et al., 2006) and the 239 signficantly, differentially 
expressed genes in Colonies 1 and 2 (FDR < 0.01). Table S11 117 genes 
that were significantly regulated by Varroa infestation (Alaux et al, 201 1) 
and immune challenge (this study). Table S12 8 genes that were 
significantly regulated by immune challenge in honey bees (this study) 
and in Drosophila (De Gregario et al. 201 1 and Roxstrom-Linquist et al. 
2004). 
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